Iterating Brownian motions, ad libitum 



Nicolas Curien and Takis Konstantopoulos 
December 2011 



Abstract 

Let Bi, B2, ■■■ be independent one-dimensional Brownian motions defined over the whole 
real line such that Bi{0) = for every i > 1. We consider the nth iterated Brownian motion 
Wn{t) = B„(B„_i(...(i?2(Si(t)))...)). Although the sequences of processes (W„)„>i do not 
converge in a functional sense, we prove that the finite-dimensional marginals converge. As 
a consequence, we deduce that the random occupation measures of W„ converge towards 
a random probability measure /ioo- We then prove that /Xqo almost surely has a continu- 
ous density which must be thought of as the local time process of the infinite iteration of 
independent Brownian motions. 

Keywords and phrases. Brownian motion. Iterated Brownian motion, Harris chain, ran- 
dom MEASURE, exchangeability, WEAK CONVERGENCE, LOCAL TIME 
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1 Introduction 

Let B+ — {B+{t),t > 0) and i?_ = {B-{t),t > 0) be independent standard one-dimensional 
Brownian motions starting from 0. The process B{t) :— B+{t) if t > and B{t) := B-{—t) 
if i < is called a two-sided Brownian motion. In this paper we study the iterations 
of independent (two-sided) Brownian motions. Formally, let Bi,B2,... be a sequence of 
independent two-sided Brownian motions and, for every n > 1 and every i € K, set 



Wn{t) := 




Burdzy [Bur9.3] studied sample path properties of the random function 1 1— > W2{t) and coined 
the terminology of (second) "iterated Brownian motion" for this object. A motivation for 
this study is that iterated Brownian motion can be used to construct a solution to the 4-th 
order PDE du/dt = ^d'^u/ dx'^; see [FuiiTtJ]. This model has triggered a lot of work, see 
[Ber96, Bur!).'^, BK!).^, BKO."), ES99, KL99] and the references therein. 

Another motivation is that the process Wn is not a semimartingale (unless n — \). Indeed, 
for n = 2, a simple calculation shows that the quadratic variation of W2 — B20B1 does not 
exist, but its quartic variation does. Similarly the 2"~variation of Wn is finite for n > 1. 
Hence, as n increases, the process Wn becomes wilder and wilder. Also, Wn is self-similar 
with index 2^", i.e., 

{Wn{at),te R) = {a^'"Wn{t),te M), 

for all a > 0. See Fig. 1 for a comparison of sample paths of Wn, n = 1,2,3. All these 
reasons make one suspect that convergence of the laws of the Wn's in a "nice" function space 
(e.g., the space of continuous functions) is impossible (see Remark 1). 
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Figure 1: First, second and third iterations of Brownian motions. 



However, following the principle of Barman saying that wild functions must have smooth 
local times (see the nice survey [GH8(I]), we prove that the occupation measures /i„ of Wn 
over [0, f ] converge in distribution towards a random measure /ioo which can be interpreted as 
the occupation measure of the infinite iteration 'Wqo". More precisely, let be the random 
probability measure defined by 



dpin / 



dtf{Wn{t)), 



(1) 



for every Borel / : M — >■ M+. Our main result is a limit theorem for the sequence (/^n)n>i- 
Restricting the integration to the unit interval is convenient and poses no loss of generality 
due to the self-similarity property of Wn- 

Let TW be the set of all positive Radon measures on E. Although we focus on the real 
line, the interested reader should consult [ , ; ■ o] for the general theory of random measures. 
Endow A4 with the topology T of vague convergence, that is, the weakest topology which 
makes the mappings 

JR 

continuous. (Here, Ck is the set of continuous functions with compact support.) A random 
measure is a random element of the space {A4,T), viewed as a measurable space with a- 
algebra generated by the sets in T. A sequence Ai,A2,... of random measures converges 
in distribution towards a random measure A if for any bounded continuous mapping F : 
{M,T) -> M we have E[F{\,)] E[F{X)] as i ^ cx3. We write 



An 



A, 



id). 



to denote this notion. Convergence of A„ to A in distribution is equivalent to: A„/ — > A/, for 
any continuous / £ Ck (see Theorem 4.2 in [KalTG]). The latter convergence is convergence 
in distribution of real-valued random variables. 

Theorem 1. There exists a random measure /loo such that /i„ -^-^ /ioo. Moreover /ioo(IR) = 
1, a.s. The random probability measure /ioo almost surely admits a density (La)aeK with re- 
spect to the Lebesgue measure such that La has compact support and is Holder continuous 
with exponent 1/2 — e for every e > 0. 

We can think of /ioo as the occupation measure of the infinite iteration of i.i.d. Brownian 
motions. Thus the random function (ia)aeR must be thought of as the local time of this infi- 
nite iteration. The convergence of the last theorem will be obtained by proving convergence 
of finite dimensional marginals of the iterated Brownian motions, namely: 

Theorem 2. For any integer p > 1, there exists a random vector {Xi)i<ci<p such that for 
any pairwise distinct nonzero real numbers xi, X2, . . . ,Xp we have the following convergence 
in distribution 

(W^"(^»))l<.<, ^ {X.)l<^<,■ (2) 
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The random variables Xi, . . . , Xp and the differences Xi — Xj, i ^ j, have all identical 
distribution which is that of a "signed exponential" with parameter 2, i.e. a distribution with 
density e"^'"^', a; G M. 

The case p = 1 has already been noticed in the Physics hterature [Tvu-04] in the context 
of infinite iteration of i.i.d. random walks. We are unfortunately unable to give an explicit 
formula for the distribution of {Xi, . . . , Xp) when p > 2, see Section 4 for a discussion and 
simulations of this intriguing probability distribution. It is quite interesting to observe that 
the marginals and differences have all the same distribution but are, of course, dependent. 
In a certain sense, the infinite iteration is both self-similar at all scales and long-range 
dependent. 

The paper is organized as follows. The first section is devoted to the proof of Theorem 2 
which is the cornerstone in the proof of Theorem 1. We then turn to the study of occupation 
measures in Section 2. The last section presents some open questions and comments. 

Acknowledgments: We are deeply indebted to Yuval Peres for insightful discussions. 

2 Finite dimensional marginals 

In this section, we prove Theorem 2. The convergence of one-dimensional marginal is a 
special case because we can explicitly give its limiting distribution, which will be of great use 
throughout this paper. In the general case, the convergence of finite dimensional marginals 
comes from ergodic property of random iterations of independent Brownian motions. 

2.1 One-dimensional marginals 

Let M* := M \ {0}. For A > 0, we denote by ±£(A) a signed exponential distribution with 
parameter A, i.e., one which has density proportional to e"'^'^', a; S M. 

Proposition 3. For any t e M* we have the following convergence in distribution 

Wn{t) ±f(2), (3) 

Remark 1. Proposition 3 already implies that the sequence of processes {Wn) is not tight 
for the topology of uniform convergence on compact intervals. Indeed, if this were the case, 
and since W„(0) = 0, we would have that sup„>]^ P{Wn(s) > rj) ^ as e ^ 0, for every 
rj > 0, and this contradicts Proposition 3. 

Proof. By standard properties of Gaussian variables we have the following chain of equalities 
in distribution 

{B,{B^{B,m)-) 
{B,{B^{^\B,{1))))..) 

{B,{\t\^/^^\BmB2{l)))..) 

n-l 

n i-^'i'""' 

1=0 

where J\fi are i.i.d. standard normal variables and ± is a independent fair random sign. It is 
then easy to see that the modulus of the right hand side of the last display actually converges 
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almost surely as n — >■ oo. Indeed the series of P(| log(|A/ip )\ > i ^) is summable and an 
application of the Borel-Cantelli lemma proves the claim. Thus, if we set 

n-l 

X := lim TT IMT' a.s., 

we have the convergence Wn(t) — ^ ±X in distribution as n — >■ oo, where ± is fair random 
sign independent of X. Notice that the limit does not depend on t as long as i 7^ 0. To 
identify the limit distribution, we note that X satisfies the following recursive distributional 
equation 

X = Vx-\Af\. (E) 

Iterating this equation and applying the same arguments as above, it is easy to see that it 
admits a unique fixed point as long as A" > almost surely. Guided by the result of [Tnr{)4], 
an easy computation shows that the exponential variable of density 2e~^^lxyo also satisfies 
(E) thus completing the proof of the proposition. □ 

2.2 General case 

The goal of this section is to prove Theorem 2 in the case p > 2. The convergence (2) will 
be achieved by applying arguments from the theory of Harris chains. The idea is to consider 
the random transformation which associates with any p points xi,. . . ,Xp G M the images 
. . . , B[xp) of these p points under a two-sided Brownian motion B and to show that 
independent applications of this map possess an ergodic property. That is, for any initial 
state (xi, . . . ,Xp), the distribution of {Wn{xi), . . . , Wn{xp)) converges weakly to a unique 
invariant probability measure. 

Let p > 2. Denote by W the set of p-uplets (xi, . . . ,Xp) of pairwise distinct nonzero 
real numbers. Note that if x = (a;i, Xp) £ TZP then its image i3(x) := {B{xi), . . . , B{xp)) 
under a two-sided Brownian motion B almost surely belongs to TV. 








Figure 2: A 6-uplet and its image after applying an independent Brownian motion. 



Proposition 4. For any p > 2, there is a unique probability measure Vp on TV such that 
if (Xi, Xp) is distributed according to Up and if B is an independent two-sided Brownian 
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motion then we have the equality in distribution 

(X,)i<,<p. (4) 
Furthermore, for any x — (xi, . . . , Xp) G TZ^ we have the following convergence in distribution 

Notice that this argument does not give an expHcit expression for the stationary proba- 
bihty measures Vp but characterizes them uniquely by (4). Let p > 2 be a fixed integer. Fix 
also a p-uplet x = (xi, . . . , Xp) £ W. For every n > 1 we set 

Wn := {Wn{xi),...,W^{Xp)). 

Thus the process (yVn)n>o is a Markov chain with state space TIP , starting from Wo := x, 
with transition probability kernel given by 

9{y- A) ^ E[{B{y^), . . . ,B(vp)) £ A], 

for any y € TV and A C TIP , where _B is a two-sided Brownian motion. We denote by i?x['] 
the expectation of this chain started from x S TIP . Since we restricted ourselves to TW ^ 
the chain (yV„) is easily seen to be irreducible with respect to the p-dimensional Lebesgue 
measure on TiP and aperiodic. We will show that the chain is in fact positive Harris recurrent, 
which will imply that it admits a unique invariant probability measure, and Proposition 4 
will directly flow from it, see [MT09, Theorem 13.0.1]. The key to prove this is to consider 
the following sets: 

Definition 1. Fix M > 1 a (large) real number. We say that the p-uplet x = (xi, ...,Xp) £ 
TIP is AT -sparse if we have 

M^^ < sup l^il < M, and min \xi — Xj\ > M^^ . 

l<i<p 1<*7^J<P 

The set of all M -sparse p-uplets in TiP is denoted by Sm- 

The basic observation is that if x = {xi, ...jXp) and y = ...Typ) are two M-sparse p- 
uplets and if i? is a two-sided Brownian motion then the two random p-uplets {B(xi), . . . , B{xp)) 
and (i3(yi), . . . , B(%jp)) are mutually absolutely continuous with Radon-Nikodym derivative 
bounded from below by a positive constant Cp^M depending on M and p only. The reason is 
that the ratio of the two densities f(xi, . . . , Xp)/f{yi, . . . , j/p) is a continuous function on the 
compact set Sm- In other words, there exists two (dependent) Brownian motions B and B 
such that B{xi) = B{yi), ...,B{xp) = B{yp) with probability at least Cp^M- Put it otherwise, 
if we flx X e Sm, then 

P(y,A) >Cp3/P(y,^), 

for all y £ Sm and all measurable A C TV. Thus Ney's minorization condition holds and, 
therefore, the set Sm is a petite set in the sense of [MT09, Chapter 5]. 

Proof of Proposition 4- In order to prove that (yV„) is positive Harris recurrent, we will show 
that, for some AT > 0, the expected return time to the petite set Sj^i, by the Markov chain 
(W,i) started from x, is bounded above by a finite constant, uniformly over all x £ Sm- 

The technical tool that we use here is the so-called drift condition, see [MTO!)], for a 
Lyapunov or "potential" function V : Tl'^ — > M+. It is required that V be unbounded on Tl^. 
The drift of such a function is defined by 



DV{x) := E^[V{Wi)]~'V{x), x G 7^P. 
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Thus, D is the generator of the Markov chain. We wiU show that there exists a Lyapunov 
function V such that, for some < a, 6 < oo, 

DV{x) < -a + bl^eSM, (5) 

for every x £ TZ^. More precisely, we show that the preceding condition is satisfied for the 
function V defined by 

V{xi,...,Xp) := max + \xi - Xj\~^^'^ , 

l<i<p — 

0<i<j<p 

and for some a,b,M > 0. In the last display, we use the notation xq := to avoid extra 
terms in the definition of V. It is convenient to consider the terms /7(x) := maxi<i<p \xi\, 
and G'(x) := J2o<i<j<p ~ comprising V, separately. For any A > 0, 

E^[U(Wi)] = E[ max \B{x,)\] = y/X E[ max \B{xi/X)\] 

l<i<p l<i<p 

Letting A = U{x.), we thus have 

E^[UiWi)]<Ciy/U{^, 

where Ci — E[max{\B{t)\ —1 < t < 1}] . For the other term, we have, with J\f a standard 
normal variable, 

i?x[G(Wi)]- J2 E[\Bix.)-Bix,)\-'/']= J2 \^.'^j\-'^'Em~'/'] 

0<i<j<p 0<i<j<p 

<C2VW), 

where C2 — pE[\J\f\~^^'^]. Putting the terms together, we find 

E^[ViWi)]<C3^V(^, 

where C3 = 2max(Ci, C2). Thus, for all x e W, 

W(x) <C3v/nx)-^(x). 

If X e Sm then y(x) < M + {p + 1)^^^ < 00. If x ^ Sm then \xi\ > M for some i 
or \xi — Xj\ < 1/M for some < i < j < P- In the first instance, V{x) > M; in the 
second instance, ^(x) > y/M. So, for M > 1, we have V{x.) > \/Af, for all x ^ Sm- Let 
M = max(16C|,l). Then for x ^ Sm, we have ^y(x) > 2C3, implying that DV{-x) < 
-C3^/V{x) < -C2MI/4. Thus (5) holds. Theorem 13.0.1 in [MT09] now applies, showing 
positive Harris recurrence. This finishes the proof of the proposition. □ 

Remark 2. The proof actually shows that there are constants C4, C5 > such that 

w(x) < -C4yy(x), 

when y(x) > C5, enabling one to establish a rate of convergence towards the stationary 
distribution. We shall not pursue this further herein. 

2.3 Properties of the z/p's 

We can think of Vp as a measure on the product space M.P (by letting it have mass out- 
side TZ^). The probability measures {i^p)p>i are consistent. This follows from the fact 
that Vp is the limit in distribution of (W„(a;i), . . . , W^„(xp)), whereas j/p+i is the limit of 
{Wn{xi), . . . ,Wn{xp),Wn{xp+i)), for any (xi, . . . , Xp, Xp+i) G TV. Therefore Vp is the pro- 
jection of Vp+i on the first p coordinates. By Kolmogorov's extension theorem, there is a 
probability measure v on M'*' such that Vp is obtained from v by projecting on the first p 
coordinates. The family {(yp)p>i,v) has some further properties. 
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2.3.1 Exchangeability 

The last statement of Proposition 4 actually shows that v is exchangeable, that is, invariant 
under permutations of finitely many coordinates. Let {Xi,X2, ■ ■ ■) be a random element of 
with law v. By the classical de Finetti/Ryll-Nardzewski/Hewitt-Savage theorem (Theo- 
rem 11.10 [Kal(j2]) we know that there exists a random element oi Ai (random probability 
measure) /ioo, such that the law of {Xi,X2, ■ ■ ■) conditionally on /j,^ is a product measure 
(the law of i.i.d. random variables). The random measure /ioo will be identified and inter- 
preted in the next section as the occupation measure of the infinitely iterated Brownian 
motion. 

2.3.2 Stationarity of the increments 

The family of measures (fp)p>i also possesses another property which can be described as 
stationarity of the increments. 

Proposition 5. Let p > 2 and let {Xi, X2, . . . , Xp) be distributed according to Vp. Then for 
every I < i < p we have 

{Xi — Xg,. . . , Xi^i — Xi, Xf+i — Xi, . . . , Xp — Xg) *■= lyp^i- 

Proof. Let {xi,X2, ■ ■ ■ ,Xp) € TZ^. It is easy to prove by induction on n > 1 and using 
elementary manipulation of the Gaussian distribution that the random vector (Wnixi) — 
Wn{xe))i^i has the same distribution as the vector {Wn{xi — xg))i^i. Notice that the vector 
{xi — xi)i^i e W"^ , thus we can apply Proposition 4 and this finishes the proof of the 
proposition. □ 

3 Occupation measure of the infinitely iterated Brown- 
ian motion 

3.1 Existence of /ioo 

Recall the definition of the random measures /i„ given by formula (1) in the Introduction, 
as well as the notion of convergence in distribution of random elements of (A^,T). With 
Theorem 2 in our hands, it is now easy to prove convergence in distribution of the random 
measures /i„. We first need a lemma, characterizing convergence in distribution of random 
elements of (A^,T), tailor-made for our case. The lemma can be of independent interest. 

Lemma 6. Let Ai, A2, . . . be random probability measures on M. The following are equivalent: 

(i) The sequence X„ converges in distribution to some random probability measure. 

(ii) For each n > 1, and conditionally on A„, let . . . be i.i.d. real-valued random 

variables each with (conditional) law A„; 

P(Xi" eAi,...,x;'eAp\ A„) = A„(Ai) • • • A„(Ap), a.s., 

for all p > 1, and Borel sets Ai, . . . , Ap. The random vector (X", . . . , ) converges 
in distribution, as n ^ 00, to some probability measure on . 

Proof. The implication (i) =^> (ii) is an easy consequence of [KaI76, Theorem 4.2]. For the 
other direction, fix / e Ck- We will show that the random variable Xnf = /r /dA„ converges 
in distribution as n 00. Consider the random probability measure 
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(i.e., the empirical distribution of {X^, . . . , X^)). We compare A„/ to d^\f ^ p-^ J2Li fi^i) 
We have ^ 

P{\Kf - U\ >e)= E[P{\Kf - U\ > e I A„)] < tt^, 

by Chebyshev's inequahty. On the other hand, fi^V converges in distribution as n — )■ oo. 
Hence A„/ converges in distribution, and thus A„ converges in distribution to some random 
measure A. It remains to prove that A almost surely has mass one. Since converges in 
distribution, it follows that the sequence of random variables {X"}n>i is tight. Thus for 
every e > there exists M > such that > Al) < e for every n > 1. Conditionally 

on A„ we have P(|X"| > M \ A„) = 1 — A„([— Af, Af]). Taking expectation we deduce that 
E[Xn{[- M , M])] > 1 - e for every n>l. This is sufficient to apply Theorem 4.9 in [Kal86] 
and deduce that A is almost surely a random probability measure. □ 

Let us go back to our setting and show that the occupation measures fj,n converge to a 
random probability measure /ioo- Let p > 1 and, conditionally on /i„, let X^, . . . ,Xp be 
i.i.d. random variables with common distribution We show that {X", . . . , X^) converges 
(unconditionally) towards the random vector (Xi, Xp) of law Vp identified in Proposition 
4. Indeed, by the definition of /l<„, for any Borel bounded function / : M^* — >■ M we have, by 
Fubini's theorem, 



i;[/(xl"\...,x("))] = E 



dui...dup /(W„(mi),..., W„(Mp)) 



dui...dup S[/(W^„(ui),...,I¥„(up))]. 



The last integral converges towards E[f{Xi, ...,Xp)] as n — >■ oo because of Theorem 2 and 
dominated convergence. Applying Lemma 6 we get the existence of a random probability 
measure fioo such that ^„ — > /ioo in distribution as n ^ oo. 



3.2 Support of /ioo 

Proposition 7. Almost surely, the random probability measure /ioo has a bounded support. 

In order to prove the last proposition we use a very general fact: 

Lemma 8. Let X„ be a sequence of random probability measures converging in distribution 
towards a random probability measure Aoo • Suppose that the support of A„ is contained in 
[An,B„\ and that the sequences of random variables (A„) and (Bn) are tight. Then Aoo has 
a compact support, almost surely. 

Proof. Let us argue by contradiction and suppose that Aoo has probability at least e > of 
having an unbounded support. By the assumption made on An and Bn there exists Af > 
such that P{\An\ > M) < e/10 and P(|B„| > Af) < e/10 for every n > 1. For the Af > 
chosen above, there exists 5 > such that we have Aoo(] — M, Af['^) > S with probability at 
least e/2. Now choose p > S^^. Suppose that, conditionally on Aoo, the random variables 
X^, . . . , X^ are i.i.d. with common law Aoo- We then have 

P ( sup \X°°\ > m] ^ E 

\l<i<p J 

> |(l-(l-<5n > £(l-e-i)/2. (6) 



PI sup |X= 

a<i<p 



> Af 



Aoo(]-M,A//n X5 
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On the other hand, if X", . . . , X" are i.i.d., conditionally on A„, we have 



P sup |Xf I >M\ < P(|A„| > M) + P(|B„| > M) < e/5. (7) 

Vl<i<p / 

Since for p fixed we have {X")i<i<p — > (^f°)i<i<p in distribution as n ^ oo, comparing (6) 
and (7) leads to a contradiction. □ 

We further establish a result on the limit of the oscillation of Wn on an interval, as 
n — oo. Recall, that the oscillation of a function / on an interval J is defined by osc(/; J) := 
suPs,tgj - /(s)|, and let, for t > 0, 

A„(t) :=osc(M^„;[0,t]). 

Lemma 9. Let D := osc{B; [0, 1]) be the oscillation of a Brownian motion on a unit interval 
and let Di, D2, . . . be i.i.d. copies of D. Then, for all t > 0, 

00 

n— >-C30 

Proof. Let 

Init) := inf Wn{x), Sn{t) := sup Wn{x). 

xelo,t] x£[o,x] 

Then A„(0 ^ 5„(t) - /„(t), and 

A„+i(i)= sup Bn+iiWn{x))^ inf Bn+iiWnix)) 

0<x<t Q<x<t 

= sup Bn+i{u) - inf B„+i(u) 

/„(t)<u<S„(t) /„(t)<u<S„(t) 

= osc{B„+i;[Init),Sn{t)]) 

^:^^osc(i3„+i;[0,A„(t)]) 

= v/A„(t) osc(B„+i;[0,l]). 

Thus iterating this equation, we get in the spirit of the proof of Proposition 3 the following 
equality in distribution 

n 

A„w'^^-"n^r'"- 

1=1 

An argument similar to the one used in the proof of Proposition 3 shows that the right-hand 
side of the last display actually converges almost surely as rt — > 00, to the infinite product 
rii^o ■ Hence A„(i) converges in distribution to the same random variable. □ 

Remark 3. Roughly speaking, the oscillation of the infinitely iterated Brownian motion is 
the same, in distribution, over any interval of any length. 

Proof of Proposition 7. By Lemma 9, A„(l) converges in distribution, it is tight. Thus the 
support of /i„ is contained in a compact interval whose endpoints are tight. Applying Lemma 
8, we deduce that, almost surely, /ioo has bounded support, as required. □ 



3.3 Density of /ioo 

This section is devoted to the analysis of the properties of the density of /ioo ■ We will proceed 
in two steps. First, using standard technique of Fourier analysis for occupation densities (see 
e.g. [Bcr69]), we will prove that /ioo almost surely has a density which is in L^. At the same 
time, we obtain some estimates about this Fourier transform. We will then use a very general 
result of Pitt [Pit 78] on local times to prove that this density is in fact continuous and even 
Holder continuous with exponent 1/2 — e for every e > 0. 
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3.3.1 Harmonic analysis of /Xqo 
For C e M, let 

*(0 := / d^oo(a;)exp(i^x), 

be the Fourier transform of the random probability measure //qo ■ 
Proposition 10. For any £^ £ R we have 

E[\mr ^ 



4 + ^2' 

Proof. By the definition of the Fourier transform <&(^), we have 



(8) 



E[\m\']=E 



By the convergences already established we have 



E 



dMoo(^)dA*oo(2/)e'«("-^) 



= lim E 



exp {l^{Wn{Ui) ~ Wn{U2))) 



where Ui,U2 are two independent random variables uniformly distributed over [0,1] and 
also independent of the sequence of Brownian motions Bi, B2, ■ . ■■ By stationarity of the 
increments (see the proof of Proposition 5), for any s,t G M, we have Wn{t) — Wn{s) = 
Wn{t ~ s) in distribution, and thus. 



E 



= hm E 



lim 

n— J-oo 



exp {imn{Ui-U2))) 



dudvE [exp(i^W„(u - v))] . 



Applying Proposition 3 and the dominated convergence theorem we get that 

4 



E 



= i?[cxp(*e[±£(2)])] - 



□ 



In particular, applying Fubini's theorem we deduce from (8) that 

2 



£;ni<i>ii^i = 



d^ < 00, 



and thus that $ e almost surely. By standard results on Fourier transforms, this implies 
that, almost surely, ^00 has density (ia)aGR with respect to the Lebesgue measure which is 
in L^. Notice that the estimates (8) in fact give us a bit more than that, namely, for every 
< s < 1/2 we have 



l<f| 



de(i + e2)1$(0P < 



a.s. 



Applying the standard Sobolev inequality, we deduce that for every < s < 1/2 we have 

||i||LV(i-2=) < C||$||//=, 

where C is a universal constant; see Theorem 1.38 in [^jC'Dll]. Since s can be arbitrarily 
close to 1/2 and since L is itself a density (thus L e a.s.) we get that L e for any 
1 < q < 00 almost surely. 
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3.3.2 Continuity of the density 

The idea to prove continuity of the density (La)aeE is to apply once again a Brownian motion 
"on top of /ioo" and to use the general theory developed in [PitTS]. Formally, if //oo has density 
(La)agR (which is in every L'', 1 < g < oo) we define the random measure /loo by 

Jloof ■■= J dtiooix) f{B{x)) = J daLa f{B{a)), 
where B is an independent two-sided Brownian motion. Clearly, we have 

Moo — i-^oo 

(equality in distribution). So it suffices to show that /loo has a continuous density. This 
follows from the lemma: 

Lemma 11. Let g : M — > he the density of a probability measure supported on a compact 
interval / C M such that (/ G for every 1 < q < oo. Let also B be a two-sided Brownian 
motion and consider the occupation measure of B with respect to the measure g(t)dt defined 
as follows 

l^f = f gmtf{B{t)), 
Jr 

for any Borel positive function f . Then the random measure /t almost surely has a density 
which is locally Holder continuous of exponent 1/2 — e for any e > 0. 

Proof. Let < 5 < 1 and choose g > 1 such that (1 + 5)q* < 2 where q* is defined by 
1/g + l/q* = 1. Applying Holder's inequality we get that 

sup^dt^^^^^ < (^y^di|g(t)r) ^%up^diK-sr(W/2^ < oo. 

We can thus apply Theorem 4 of [ ] and get that /t a continuous density L over M which 
satisfies a local Holder condition \L{x) — L{y)\ < K\x — , for any 5' < 6/2. □ 

4 Comments 

Occupation over other measures. In this work we considered the occupations mea- 
sures /i„ of Wn over the time interval [0, 1] but the proofs apply in the more general setting 
of occupation measures of Wn over any probability measure on K which is not atomic. More 
precisely if A is a probability measure on M which has no atom, we define the occupation of 
Wn over A by 

mI^'/ / Xidt)fiWnit)), 

Jr 

for any Borel / ; M — > Then the random probability measures /il^^ converge towards 
/too in distribution as well. 

Exphcit finite dimensional marginals. The consistent family distributions > 
1} introduced in Proposition 4 are the limiting distributions of the finite-dimensional marginals 
of the W„'s. They are characterized by (4). Although they arise naturally in the study of 
iteration of Brownian motions, to the best of our knowledge, they have not been investigated 
so far for p > 2 and in particular no explicit formulae are known for the density oi Vp, p > 2. 

In Fig. 3 many independent points have been sampled according to 1/2- A clear shape "sea 
star" emerges and we conjecture that the level-lines of the density of V2 would be dilatation 
of this unique shape. However we do not have any candidate for this density. 
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Figure 3: About 300000 independent samples of the distri- 
bution 1/2- Simulation realized by Alexander Holroyd. 

Funaki's PDE. Following Funaki [Fun79], we can show that Wn solves the PDF du/dt = 
d^u/dx^, where k = 2", with initial data u{t = 0,x) = uo{x), x & R, where uq is a func- 
tion satisfying growth conditions [Fun79]. The solution is given by Un{t,x) — E[uo{x + 
Wn{2'^ ~^t))]. Although it is probably not possible to give a PDE meaning to the exchange- 
able collection of random variables {Wooix))x,£R, we remark that the limit to the signed 
exponential translates into lim„_j.oo w„(t/2^ = J-^Uo{x — y)e~'^^^^ dy . 

Ray-Knight theorem. In the spirit of the famous Ray-Knight theorems, do we have 
another way to describe the density of /ioo as a diffusion process? 

Reflected Brownian motion. It is possible to extend part of the previous work to 
iterations of reflected Brownian motions. Namely, it is likely that Proposition 4 goes through 
but Proposition 5 fails, and thus the analysis of the Fourier transform of the limiting measure 
would be more challenging. 

Fractional Brownian motion. Fractional Brownian motion generalizes Brownian mo- 
tion in that it is a Gaussian i?-self-similar process with stationary increments, where < 
H < 1. Again, it is likely that Proposition 4 go through, but the limit is even harder to 
describe. Of course, this raises the following interesting question: what kind of processes can 
be iterated ad libitum and result into some kind of limit? 
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